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Abstract 

We derive and employ a semiclassical spherical flow model to describe the 
atomic dynamics and calculate the current-induced force in the presence of electrical 
current. The conduction electrons serve as an ideal gas in classical mechanics flowing 
over atoms while using quantum mechanics for calculating the fluid pressure and 
viscosity. Two new kinds of current-induced force are derived based on present model, 
which are electron density force and electron viscous force. The calculated results are 
in good agreement with previous experimental results [Science 328, 736 (2010)]. Our 
model is also suitable for the motion and force analysis of nanoislands and voids. 
Combined with previous results, we firstly build the connection between 
electromigration and electroplasticity and help solve some long-standing puzzles. 

Introduction. — Electroplastic effect was discovered by Troitskii and Likhtman 1 . 
It was found that there was an immediate stress drop associated with application of a 

1 2 

current during plastic deformation. ' The stress drop during electric current 
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application was termed the electro-plastic (EP) effect. This indicated that electric 
current might have a positive effect on material plasticity during plastic deformation, 
which had been broadly studied during the last two decades. 3 " 19 In particular, there 
was an optimum strength of electric field which decreased the flow stress 
approximately 20%. 4-5 and the mobility of metal defects (including dislocations, voids, 
grain boundaries) was found increased due to the applied large current. 3 ' 6 ' 7 ' 16 " 19 
However, there were still some puzzles, such the mobility of dislocation increased 
both with or against the current direction during current application, 16 " 19 which could 
not be well explained. On the microscopic scale, the large current in a metal can 
generate preferential diffusion of atoms in one direction. This phenomenon was called 
electro migration. 20 ' 21 The effect of current-induced forces in atomic and molecular 

22 27 

wires has been studied in recent work. " Though the current-induced force had been 
broadly measured or calculated, however, there's still a huge difference between the 
predicted value and the experimental results due to the complexity of many body 
problems. Also, some phenomenon cannot be well explained, such as the grain 
boundary drag effect, 28 the motion of nanoislands. 29 A connection between them is 
also fail to be built due to topological problem (dislocation) and the complexity of 
atomic structure. 

A unified model including all aforementioned effects is highly desirable for 
further exploration in this field. In this article, the conduction electrons are degenerate, 
free, no n- interacting particles, which serve as an ideal gas in classical mechanics, with 
the pressure and viscosity determined by quantum mechanical effects. Two new kinds 
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of current-induced force are derived and they can be roughly represented as 3p e haV 
and Ap e K . Despite its seemingly simplicity, the calculated results agrees surprisingly 
well with previous experimental results, which is even much more accurate than DFT 
calculation results. These two kinds of force can also be applied to different models to 
explain some puzzles. For example, combined with Frenkel-Kontorova model, the 
connection between electromigration and electroplasticity is firstly built. The 
dislocations are driven by the force that is exerted on atoms, which is completely 
different from conventional deformation mechanisms and help explain some puzzles. 

Electron density force (EDF). — The atom movement in the electron gas is 
reduced to the classical movement of a sphere in unbounded ideal fluid. We mainly 
analyze the movement of atoms in dislocation line. The atom in the dislocation line 
(hereafter called dislocation atom) is the sphere, and the electron gas is the unbounded 
ideal fluid. The fluid is the relative reference frame, and the sphere moves in the fluid 
with velocity V (V = v ato -v e ) , where v ato is the atom velocity, v e is the free 
electron drift velocity. On the sphere surface, the distribution of relative velocity is 
v r=a = — V sin 6 , where 0 is the angle between the radius and x-coordinate. 30 
Average the horizontal component of velocity of the left side of sphere, we arrive at 



Similarly, the velocity of the other side can be got as —V. Assuming that A is at 
the left side of dislocation and B is at infinity. The mean horizontal velocity of the 
electron gas at A is v A , the horizontal velocity of the electron gas at B is v e , which is 




(1) 
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the electron gas drift velocity. The pressure of the electron gas at A is P A . The 
pressure of the electron gas at B is P B = P and p B = p , where P is the mean 
electron gas pressure, p is the mean electron gas mass density. Each side of the 
dislocation atom is assumed to meet Bernoulli equation with different density, 
respectively. The pressure equation for electron gas at A is obtained as: 



2 2 



P A =^P + p A Ve J A (2) 
p 2 



Similarly, for the other side C 



2 ..2 



P c= ^ P+Pc ^ (3) 
p 2 

3 3 1 

where v c =v A = — V + v alo =—v e +—v at0 . According to D'Alembert paradox, the 
pressure of both sides should be the same. However, the electron gradient between 
two sides of the atoms in the dislocations line leads to the electron density force. The 
direction of EDF is only determined by the electron density gradient. Then, the 
expression for EDF can be obtained as 



2 

V. 



3 1 

— V H V 

4 e 4 a >° 



°ed F =P A -P c =^P + Ap ^— '- (4) 

p 2 

where p = p e -m e , p e is the electron distribution density, m e is the electron 
mass. Then the basic equation for EDF is obtained as follow 



2 



r 3 1 - : 

— V H V 

4 e 4 0,0 



P + A Pe m e ^— (5) 

Pe 2 

For conduction electrons, only those very close to the Fermi energy can 
participate in the conduction process. Therefore, their pressure can be got by small 
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deduction as P = —p e E f where E f is Fermi energy level. For the second part of 
equation (5), it derives from the kinetic energy from electrons to ions. It can be 
replaced by 

/ 

P e v e v sca Tcoi E f = P e v e cr sca -^E f =E f (6) 

v 

e 

where a sca is the scattering cross-section of ion, r col is the mean time interval 

between two successive collisions, l mfp = ( p e cr sca ) ' = p~ m is the mean free path 

between electrons and ions. Assume the atom velocity in current direction is zero. 

Then equation (5) can be rewritten as 

a EDF = -Ap e E f +^~A Pe E f = OMAp e E f (7) 
5 16 

Ap e is the electron distribution density gradient between two sides of 
dislocations. Electron density force is like the nature wind which is created by the air 
density difference/gradient. The distribution of electron is significantly affected by 
various factors such as charged vacancies and congregated solute atoms. 31 " 35 Between 
the two sides of dislocation, under the effects of charged vacancies, a vivid contrast of 
electron distribution density is formed. The electron density gradient has been 
observed in previous researches. The constrictions in current flow path can increase 
local density up to about a factor of 2. 36 Also, atomic relaxation at step edges 

37 38 

redistributed local charge density, ' with enhancement up to a factor of 10 on 
Cu(532) near kink sites. 39 

An equation for roughly estimating dislocation charge polarization is proposed 
based on the assumption that gliding dislocations sweep up charged vacancies. 40 For a 
dislocation pileup of 2 / 0 length locked at x = ±/ 0 , the charge distribution density can 

5 
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be written as follows 



Pchar(x,t)=A 0 Qxp(-rt 2 




x 2 +a 2 (t) (x + lj +a 2 (t) (x-lj +a 2 (t) 



2 1 1 



(8) 



where A 0 is the amplitude of the vacancy charge distribution at initial moment 
( t = 0 ); r is supposed to be a constant. The dislocation is at initial moment t = 0 , 
the density equation is thus 



The dislocation pileup consists n dislocations, and the dislocation is equidistant 
from each other. Then, the x coordinates of the top dislocation and the next 
dislocations are x l = /„ and x 2 = / 0 - 2/ 0 In , respectively The charge density 
gradient Ap char , which is proportional to Ap e , can be obtained as 



Equation (10) builds a connection between plastic deformation and electron density 
gradient. During plastic deformation, the electron density gradient is mainly 
determined by charged vacancies or solute atoms which are swept up by moving 
dislocations. The electric field or current can change direction of electron density 
gradient, thus forming new motive force for atom or dislocation motion. This 
prediction is recently demonstrated by reference 41. Using DFT calculation, the 
charge density gradient of Ag can be represented as: 



27 3 

And the simulation results indicates that could be of 4X10" m" on chain and 
dislocation lines, which is directly suitable for present model. 

6 




(9) 



Ap e cc A Pchar =p(l 0 ,0)- p(l 0 -2l Q I n, 0) 



(10) 




(11) 
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The current-induced stress is about a few GPa. 7-913 In reference 41, Ap e is 
firstly calculated using first principles recently and Ap e could be of 4xl0 27 m 3 . 

28 3 

The conduction electron density p e for single-crystal Ag is 5.85 x 10 m" , Fermi 
energy is 5.48eV, therefore, according to equation (7), <j edf = 2.94Gpa . Equation (7) 
contains two parts. The left part derives from the electron density gradient, while the 
other part derives from the quantum kinetic energy transfer. We use classical fluid 
mechanics to deduce the basic equations for EDF. Then, we further transform it to 
quantum mechanical forms, and the calculated results demonstrate the success for the 
transformation we have made. 

Electron viscous force (EVF). — The other kind of current-induced force is caused 
by electron viscosity. In classical fluid mechanics, the viscous drag formula has been 
given by Navier-Stokes 30 

F = 6xpaV (12) 
where V is the relative velocity of the sphere to the fluid and V = v ato -v e , u 
is dynamic viscosity, a is the sphere radius. We assume the conduction electrons as 
degenerate matter (Fermi gas). A Fermi gas of weakly coupled systems obey 
universal hydrodynamics. 42 ' 43 The viscosity has a universal quantum scale p e ir , where 
it is Planck's constant h divided by 2n . A more accurate form of this viscosity 
can be got as 

M = <*„Pe+ ( 13 ) 

where a p is a universal dimensionless coefficient, the electron viscous force can 
thus be rewritten as 

7 



Generated by Foxit PDF Creator © Foxit Software 
http://www.foxitsoftware.com For evaluation only. 

3a M p e haV (14) 
It is worth noting that equation (13) is applicable to conduction electrons when 
the mean free path / - of electrons is much smaller than the characteristic size R 
of the system. The so-called Knudsen number is K = I. I R , 42 For conduction 

electrons, the characteristic size is determined by the de Broglie wavelength 

A = h/^2m e E f and l emfp = p~ m . /A = 0.49. 

29 

We then choose the value in reference to make a comparison between the 

predicted EVF value and the measured one. According to v e = , when j is 

Pe e 

9 2 

6.7><10 A/m , v e is 0.715 m/s. The ratio of mean energy per conduction electron to 
Fermi energy E I E f is about 1, according to reference 42 , the viscosity coefficient 
a is 2.30-2.50, then, the dynamic viscosity is /u = (1.42 ~ \ .54)x\0 5 pa -s . The 
atom radius of Ag + is 1.26xl0" 10 m, then F EVF Ij = (2.25 ~ 2.44) xl0 _3 e Fern/ A . 
Compare with the measured value F I j = (2.2 + 0.9) x 10 ^ eVcm I A , 29 However, the 
expected values based on previous models using first principles are 
3.0x10 ^ eVcm I A and 6.6x10 5 eVcm/ A , 35 ' 36 and there's still lack of an accurate 
and effective method for calculating the current-induced force at present 29 We 
introduce fermions viscosity into the classical stokes drag formula and it's surprising 
that the calculated value is of such high accuracy. 

Chenggang Tao et al visualized the electron scattering force in Nanostructure by 
observing the current-biased displacement of monatomic islands of several 
nanometers on single-crystal Ag. 29 It was found that the velocity of the islands varies 
inversely with radius, which is contradictory to the prediction of pure 
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attachment/detachment mechanism. However, this phenomenon can be well explained 

by our theorem. Assuming the nanoisland is a disk, then the drag force of the disk is 
32 

F disk = - — juVR 0 , 30 the island area density is p is , the island radius is R 0 , the island 



drift velocity is v 0 . 



dv 0 = 32/yfr(v g -v„) 
dt 3p is 7rR 0 



(15) 



Then, we can get 



32/yfr f 

v e -Ae 3p ^ =v 0 (16) 



Because v 0 =0 when t = 0 , then A = v e , therefore 

32p c * 



v 0 =v e (\-e 3 ^) (17) 



dv 0 32p e -hv e - 



(18) 



dR 0 3p is x% 

According to equation (18), the island drift velocity v 0 varies inversely with 
radius R 0 , which well suits the experimental results. Additionally, according to 
equation (18), the island drift velocity will never be larger than the electron drift 
velocity. This phenomenon is thus explained in a completely different way from the 

28 

diffusion mechanism. An important observation was made by K.C. Chen et al. They 
found that the grain boundary slowed down electromigration rate by one order of 
magnitude. According to equation (11) and (7), this time the electromigration 
direction might be opposite to the electron density gradient in grain boundary, which 
leads to the drag effect. 

Application. — While previous investigation had proposed the concept of 'electron 
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wind force' which might help promote the deformation, 4 ' 10 they failed to demonstrate 
that electric current could exert force on topological defects (dislocation). In 
Frenkel-Kontorova model, the dislocation motion which is driven by dc (direct 
current) force can be represented as 44 

iij (t) + rju, (t) + smu, ~g(u M + u,_ i -2m,) = F (19) 
where the dot stands for the time derivative, 77 is the damping coefficient, u, is the 
coordinate of the / - th atom, g is the elastic constant, F is the applied dc force to 
all atoms. The dc-driven FK model was investigated by a series of papers by Braun et 
al. 45 " 47 A similar ac-driven dislocation motion has also been given by previous 
researches. 48 " 49 However, while the FK model demonstrate the effect of dc current on 
dislocations, it is lack of an effective representation of dc force. Combined with our 
model, dc force F in equation (19) can be replaced by equation (14) as 

ii, (t) + Tju, (?) + sinw, -g(u l+l +u,_ 1 -2u,) =3a fJ p e ha(V e -u,(t)) (20) 
according to equation (14), the plastic deformation is promoted by EVF. If the 
drive force is EDF, according to equation (7) and (11), EDF can help explain the 
increased mobility of grain boundary. 7 It had been pointed out that current-driven 
motion of the dislocation is possible if the drive's amplitude exceeds a threshold 
value. 48 ' 49 This, combined with equation (20), can help explain why there is an 
optimal strength of electric current 4 " 5 and can help define the threshold value of 
current density. The electron wind force-driven dislocation motion had been directly 
observed and there existed an threshold value of the applied electric field for driving 
these dislocations. 50 According to equation (20), EVF and EDF provide extra motive 

10 
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force for dislocation motion. They thus help depin dislocations and reduce 
double-pinned Frank-Read sources, thereby reducing dislocation density. This 
explains the dislocation density reduction effect. 14 ' 15 

In reference 51,52, the solution has been obtained for the motion of void driven 
by electric potential and gradient stress field. If we consider the effect of EDF, which 
results from the electron density gradient of the void, the chemical potential of the 
atoms on the interface should be 

M = H> -Cly s K-Z*eU-Cla mF -Qa„ +QU S (21) 

where ju 0 is the reference value of the potential, y s is the interface energy, k 

is the curvature of the interface, Q is the atomic volume, U is the applied electric 
potential, U s is the strain energy density of matrix material, <j n is the normal stress 

acting on the interface. The atomic flux along the interface is given by 



J = - 



DAdju DA 



( 



7s 6 ^-?eE i + d ^ + d ^-™\ (22) 

OS OS OS OS J 



QkT ds KT 

where Z> is the interface diffusivity, 8 j is the thickness of the interface, kT 
has the usual meaning, s is the arc length of the interface, E j =-V \U . The mass 
conservation requires that dJ I dy = V / Q , then integrate it from y to 



y = d^{\-m)l (l + m) , where m is the shape factor, we get 

f n ~\ 



dJ—-y 
Vl + m 



+ C (23) 



where the constant C can be determined by boundary condition, m is the 
shape factor, then 



C = 



kT 



V OS OS J 
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(24) 
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where — depends on electron density gradient. Under the effect of 

ds 

gradient electron density gradient, atoms diffuse along the interface from one side of 

high density to another of lower one. Therefore, — = -V s cr ££1F < 0 .Symmetry 

ds 

requires that the flux at y = 0 must be zero, thus 

, = 5M^U +V „-^ + ^) ,25) 
dkT V 1 - m \ ds ds ) 

According to equation (25), electron density force contributes to the motion of 
voids. 3 ' 6 Additionally, EDF helps explain the variation of the void velocity which 
might result from the variation of electron density distribution. 3 ' 6 It had been shown 
that dislocations increase their velocities when moving both with and against the 
current and the dislocations move along current direction approximately 10-20% 
faster than against it, which is called non-polarization effect. 16 " 19 The improvement of 
dislocation mobility is the result of various effects including EDF and EVF and the 
other part in equation (25), such as the gradient stress field. 51 ' 52 the plasticity can be 
improved even when EVF is adverse. 

Conclusion. — The current-induced force had been called "electron wind force" 
for a long time, however, it was merely a concept. Now, in present model, this force is 
finally demonstrated as a real electron wind. Our model assumes the electron gas as 
Newton fluid, which is a very bold assumption and completely different from 
previous researches. However, the calculated results and well explanation of relative 
phenomenon strongly corroborates our assumption. We calculate the electron force 
using Fermi viscosity, Fermi pressure, Navier-Stokes formula and Bernoulli formula, 
and the calculated value agrees well with experimental value. We have also firstly 

12 
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built the connection between electromigration and electroplasticity combined with 
previous results, by which some long-standing puzzles are solved. Our model is also 
suitable for the motion and force analysis of materials of various scales, such as atoms, 
nanoislands, dislocations and voids. 
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